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Objectives:  This  ultimate  aim  of  our  work  is  to  develop  new  concepts,  theory,  and  algorithms  for  control 
and  signal  processing  using  ideas  inspired  by  molecular  systems  biology.  Cell  biology  provides  a  wide 
repertoire  of  systems  that  are  strongly  fault-tolerant,  nonlinear,  feedback-rich,  and  truly  hybrid,  while  mak¬ 
ing  effective  use  of  highly  heterogeneous  sensing  and  actuation  channels.  Advances  in  genomic/proteomics 
and  molecular  systems  biology  research  are  continually  adding  detailed  knowledge  of  such  systems’  archi¬ 
tecture  and  operation,  thus  offering,  in  principle,  a  powerful  source  of  inspiration  for  innovative  solutions  to 
problems  of  control  and  communication,  sensor  and  actuator  design,  and  systems  integration. 

Status/Progress: 

During  the  final  year  of  this  grant,  the  PI  and  his  collaborators  largely  completed  the  research  started  in  the 
first  year,  resulting  in  significant  new  results,  which  were  documented  in  submitted  papers,  on  the  formu¬ 
lation  of  an  approach  to  stochastic  monotone  control  systems  and  on  the  development  and  application  of 
tools  including  contraction  theory  to  the  study  of  stability.  The  effort  during  the  grant  period  involved  the 
PI  as  well  as  graduate  students  Zahra  Aminzare,  Maja  Skataric,  and  Michael  de  Ereitas,  the  latter  two  of 
whom  completed  their  Ph.D.  theses  in  2014.  The  thesis  by  de  Ereitas  completed  the  work  started  early  on 
in  this  project  and  provided  a  general  small-gain  theorem  for  random  dynamical  systems  with  inputs  and 
outputs,  motivated  by  bio-molecular  systems.  Work  with  Aminzare  provided  new  results  on  contractions 
for  distributed  and  interconnected  systems,  and  was  motivated  by  a  molecular  binding  model.  Work  with 
Margaliot  and  Tuller  developed  new  notions  of  “weak”  contractions,  motivated  by  problems  from  protein 
translation  machinery.  Work  with  Skataric  and  Nikolaev  dealt  with  the  use  of  singular  perturbation  tools  to 
show  that  scale  invariance  is  subject  to  fundamental  limitations,  again  motivated  by  a  set  of  questions  arising 
from  a  concrete  biological  context. 

A  brief  technical  discussion  of  recent  work  follows.  The  previous  reports  for  years  1  and  2  should  be 
consulted  for  discussion  of  work  conducted  previous  to  the  start  of  the  last  period. 

Contraction  theory 

As  described  in  previous  reports,  contraction  theory  is  a  powerful  tool  for  analyzing  nonlinear  dynamical 
systems.  A  system  is  contracting  if  all  the  trajectories  in  some  region  of  the  state  space  converge  to  each 
other  at  an  exponential  rate.  This  implies  many  desirable  properties  including  convergence  to  an  attractor, 
and  entrainment  to  a  periodic  excitation.  Contraction  theory  has  found  numerous  applications  in  control 
theory,  observer  design,  synchronization  of  coupled  oscillators,  and  more.  Recent  extensions  include  no¬ 
tions  of  partial  contraction,  analyzing  a  network  of  interacting  contracting  elements,  LaSalle-type  principles 
for  contracting  systems,  Lyapunov-like  characterization  of  incremental  stability,  and  applications  of  pseudo¬ 
contraction  to  study  contraction  to  a  subset  of  the  state-space  (e.g.,  the  set  of  states  with  equal  coordinates). 
We  have  continued  to  study  many  of  these  topics. 

Two  papers  have  been  published  or  submitted  during  this  period:  [8]  and  [6],  in  which  we  studied  a  notion 
of  “weak”  contraction  that  allows  contraction  up  to  an  arbitrarily  small  transient  in  time,  and  applied  these 
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ideas  to  the  entrainment  of  translation  elongation  to  intra-eellular  signals  sueh  as  tRNAs  levels  and  other 
faetors  affeeting  translation. 

Motivated  by  an  important  bioehemieal  example,  we  developed  new  theorems  regarding  the  synehronization 
of  agents,  intereonneeted  by  diffusive  links,  in  a  journal  paper  [3]  under  review  (a  eonferenee  version,  with 
a  summary  as  well  as  additional  material,  will  appear  in  [4]).  In  these  papers,  we  diseuss  the  applieation  of 
eontraetion  to  synehronization  of  diffusively  intereonneeted  eomponents  deseribed  by  nonlinear  differential 
equations.  We  provide  estimates  of  eonvergenee  of  the  differenee  in  states  between  eomponents,  in  the 
eases  of  line,  eomplete,  and  star  graphs,  and  Cartesian  produets  of  sueh  graphs.  We  base  our  approaeh 
on  eontraetion  theory,  using  matrix  measures  derived  from  norms  that  are  not  indueed  by  inner  produets. 
Sueh  norms  are  the  most  appropriate  in  many  applieations,  but  proofs  eannot  rely  upon  Lyapunov-like  linear 
matrix  inequalities,  and  different  teehniques,  sueh  as  the  use  of  the  Perron-Frobenious  Theorem  in  the  eases 
of  or  L°°  norms,  must  be  introdueed. 

In  [1]  we  studied  the  speeial  ease  of  norms  and  generalizations  to  reaetion-diffusion  partial  differential 
equations.  We  also  prepared  a  long  expository  paper  on  eontraetion  theory  and  applieations  to  nonlinear 
stability  for  presentation  at  an  invited  Tutorial  at  the  upeoming  2014  IEEE  Conferenee  on  Deeision  and 
Control  [2].  This  paper  provides  a  self-eontained  introduetion  to  some  of  the  basie  eoneepts  and  results 
in  eontraetion  theory,  diseusses  applieations  to  synehronization  and  to  reaetion-diffusion  partial  differential 
equations,  and  poses  several  open  questions. 

Random  dynamical  systems 

In  the  1980s,  E.  Arnold  eoneeived  an  elegant  and  deep  approaeh  to  the  foundations  of  random  dynamies.  His 
paradigm  of  a  Random  Dynamical  System  (RDS  for  short)  is  based  on  an  ultimately  simple  idea:  view  an 
RDS  as  eonsisting  of  two  ingredients,  a  stoehastie  but  autonomous  “noise  proeess”,  and  a  elassieal  dynami- 
eal  system  that  is  driven  by  this  proeess.  The  noise  proeess  is  deseribed  by  a  measure-preserving  dynamieal 
system.  It  is  typieally  probabilistie,  representing  for  example  environmental  perturbations,  internal  variabil¬ 
ity,  randomly  lluetuating  parameters,  model  uneertainty,  or  measurement  errors.  But  the  formalism  allows 
for  deterministie  periodie  or  almost-periodie  driving  proeesses  as  well.  The  resulting  theory,  developed 
sinee  by  many  authors,  provides  a  seamless  integration  of  elassieal  ergodie  theory  with  modern  dynamieal 
systems,  giving  a  theoretieal  framework  parallel  to  elassieal  smooth  and  topologieal  dynamies  (stability, 
attraetors,  bifureation  theory,  and  so  forth),  while  allowing  one  to  treat  in  a  unified  way  the  most  impor¬ 
tant  elasses  of  dynamieal  systems  with  randomness  — ^random  differential  or  differenee  equations  (basieally, 
deterministie  systems  with  randomly  ehanging  parameters),  or  stoehastie  ordinary  and  partial  differential 
equations  (white  noise  or,  more  generally,  martingale-driven  systems  as  studied  in  the  ltd  ealeulus).  In  pre¬ 
viously  reported  work,  we  introdueed  a  notion  of  random  dynamieal  systems  with  inputs,  providing  several 
basie  definitions  and  results  on  equilibria  and  eonvergenee.  The  main  goal  has  been  all  along  to  propose  a 
new  RDS-based  formalism  for  random  eontrol  systems,  that  is,  systems  with  inputs  (and  outputs),  whieh  we 
abbreviate  RDSI  (or  RDSIO). 

Our  motivation  for  studying  RDS  with  inputs  and  outputs  arises  from  the  need  to  provide  foundations  for  a 
eonstruetive  theory  of  intereonneetions  and  feedbaek  for  stoehastie  systems,  one  that  will  eventually  general¬ 
ize  sueeessful  and  widely  applied  deterministie  approaehes  to  the  analysis  and  design  of  dynamie  networks. 
Intereonneetion  theory  has  proved  very  useful  for  the  analysis  of  even  deterministie  systems  with  no  inputs. 
It  has  been  very  sueeessful  when  eombined  with  tools  from  passivity  theory,  input-to-state  (ISS)  stability, 
and  monotone  systems  with  inputs  and  outputs.  The  key  need  is  to  have  an  appropriate  elass  of  systems 
that  are  CICS  (eonverging-input  eonverging-state)  eoupled  with  what  are  generieally  ealled  “small-gain  the¬ 
orems”  (essentially,  asking  that  the  feedbaek  loop  results  in  a  eontraetion  in  an  appropriate  sense).  CICS 
theorems  were  proved  in  the  previous  year  and  have  been  already  published.  The  final  sfep  in  Ibis  projeef 
was  aehieved  wifh  fhe  sueeessful  sfafemenf  and  proof  of  a  small-gain  fheorem  for  monofone  RDSI,  and  Ibis 
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was  the  main  result  of  de  Freitas’  Ph.D.  thesis,  and  is  reported  in  the  reeently  submitted  paper  [5]. 
Scale-Invariance 

We  have  largely  eompleted  our  program  of  theoretieal  developments  eoneerning  the  seale-invarianee  prop¬ 
erty  in  enzymatie  networks.  The  phenomenon  of  fold-ehange  deteetion,  or  seale-invarianee,  is  exhibited  by 
a  variety  of  sensory  systems,  in  both  baeterial  and  eukaryotie  signaling  pathways.  A  eommon  theme  in  the 
systems  biology  literature  is  that  eertain  systems  whose  output  variables  respond  at  a  faster  time  seale  than 
internal  eomponents  give  rise  to  an  approximate  seale-invariant  behavior,  allowing  approximate  fold-ehange 
deteetion  in  stimuli.  Indeed,  in  our  own  results  reported  in  previous  grant  periods,  and  already  published, 
we  argued  that  all  3-node  scale-invariant  enzymatic  networks  are  based  upon  this  mechanism.  The  paper  [7] 
established  a  fundamental  limitation  of  such  a  mechanism,  showing  that  there  is  a  minimal  fold-change  de¬ 
tection  error  that  cannot  be  overcome,  no  matter  how  large  the  separation  of  time  scales  is.  To  illustrate  this 
theoretically  predicted  limitation,  we  discussed  two  common  biomolecular  network  motifs,  an  incoherent 
feedforward  loop  and  a  feedback  system,  as  well  as  a  published  model  of  the  chemotaxis  signaling  pathway 
of  Dictyostelium  discoideum. 
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